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\^ . We have carried out mapping observations of the molecular core associated with 

the young Class protostar, IRAM 04191+1522, in the CH3OH (J K =2 K -1 K ) and C 34 S 
(J =2-1) lines using the 45 m telescope at Nobeyama Radio Observatory. Our ob- 
servations have revealed that there is a condensation associated with the protostar, 
elongated in the east-west direction mostly perpendicular to the axis of the associated 
CO outflow. Its size and mass are estimated to be 0.07 pc x 0.04 pc and 2.3 M , re- 
spectively, from the CH3OH data. In addition to the elongated envelope, two compact 
(~ 0.03 pc) condensations were found in the CH3OH line at the southern edge of the 
elongated envelope, where the blueshifted CO outflow emerging from the protostar is 
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located. In contrast to the elongated envelope, those compact CH3OH condensations 
show much larger line width (up to 2.0 km s _1 ) with centroid velocities blueshifted by 
~ 0.8 km s . The compact condensations have momenta (~ 0.06 M km s~ x ) compa- 
rable to that of the blueshifted molecular outflow. In addition, they are gravitationally 
unbound, and most probably will dissipate eventually. These results suggest that the 
compact condensations are probably formed in the course of interaction between the 
outflow and the ambient gas surrounding the protostar, and that such interaction may 
cause dissipation of a part of the ambient gas. No drastic, localized enhancement of the 
CH3OH abundance is, however, observed toward the compact condensations, implying 
that there seems to be no significant shock heating at the compact condensations in 
spite of the interaction with the outflow. This may be because the CO outflow velocity 
(< 10 km s _1 ) is too low to cause effective heating to release CH3OH on dust grains 
into gas phase. 

Subject headings: ISM: abundances — ISM: individual (IRAM 04191+1522) — ISM: 
molecules — stars: outflow — stars: formation 
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1. Introduction 

It is widely accepted that molecular outflows take place during low-mass protostellar forma- 
tion. The origin of molecular outflows is still debatable, whereas recent observational and theoretical 
studies have suggested that high-velocity neutral jets from the vicinity of protostars drive molecu- 
lar outflows by entraining the material along the jet path (Raga & Cabrit 1993; Dutrey et al. 1997; 
Gueth et al. 1998; Hatchell et al. 1999a; Hatchell et al. 1999b). Those jets could cause shocks and 
alter the physical and chemical conditions of the entrained material drastically. In fact, there is 
observational evidence of such shocks from some Class sources, such as L1157 (Mikami et al. 1992; 
Umemoto et al. 1992; Bachiller et al. 1995; Avery & Chiao 1996; Bachiller & Gutierrez 1997; Gueth et al. 1998; 
Umemoto et al. 1999; Bachiller et al. 2001), NGC 1333 IRAS 2 (Sandell et al. 1994; Bachiller et al. 1998), 
IRAS 4 (Blake et al. 1995), IRAS 16293-2422 (Blake et al. 1994; Ceccarelli et al. 2000; Hirano et al. 2001), 
L1448 (Bachiller et al. 1991; Dutrey et al. 1997), HH25MMS (Gibb & Davis 1998), and BHR71 
(Bourke et al. 1997; Garay et al. 1998). These sources show that, in the positions where jets hit 
and entrain the ambient dense gas, molecular lines such as SiO, SO, CS, and CH3OH are much 
more intense than those in the positions of the protostars. These results have been interpreted 
as the case where molecules on dust grains evaporate into gas phase through the C-type shock 
heating or sputtering due to the interaction with jets, resulting in enormous abundance enhance- 
ment of these molecules (Kaufman & Neufeld 1996a; Kaufman & Neufeld 1996b; Hollenbach 1997; 
Pineau des Forets et al. 1997; Schilke et al. 1997; Bergin et al. 1998; Buckle & Fuller 2002). 

However, studies of the interaction between outflows and parent cloud cores have been so far 
limited to some representative sources with strong shocks like those mentioned above. Hence, we 
do not know how ubiquitous and how important the interaction between outflows and cores is in 
the process of low-mass star formation. For example, there is no thorough observational studies 
on this issue in the Taurus molecular cloud complex, even though it is one of the nearest (140 pc ; 
Elias 1978) and the most representative low-mass star-forming regions. 

IRAM 04191+1522 (hereafter referred to as IRAM 04191) is a Class protostar, recently dis- 
covered by Andre et al. (1999a) in the southern part of the Taurus molecular cloud. This source, 
with a very low bolometric temperature (T^) of ~ 18 K, is located at ~ 1' southwest of the 
Class I infrared source IRAS 04191+1523 (not to be confused with IRAM 04191), and is consid- 
ered to be one of the youngest protostars among the known Class objects (Andre et al. 1999a; 
Andre et al. 1999b). IRAM 04191 is associated with a well-collimated molecular outflow that ex- 
tends from northeast (red) to southwest (blue) (Andre et al. 1999a; Lee et al. 2002). Since such 
an extremely young protostar should be surrounded by a plenty of material which is less disrupted 
by the associated outflow, we may expect a more clear case of the interaction between the outflow 
and the surrounding material. In fact, Belloche et al. (2002) have recently found evidence of such 
an interaction toward the blue lobe of the molecular outflow in the CS(2-1) line. It is therefore 
quite interesting to study the interaction between the outflow and the core toward this very young 
protostar in detail. 
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We adopted a band of CH3OH ( Jk=2k _ 1k) lines (Friberg et al. 1988) as a probe to study inter- 
action between the outflow and the parent cloud core associated with IRAM 04191. As mentioned 
above, the CH3OH line can be a good tracer of shock interaction, since CH3OH molecules on dust 
grains evaporate into gas phase by mean of shock heating and the abundance is greatly enhanced to- 
ward those interaction regions (Sandell et al. 1994; Bachiller et al. 1995; Blake et al. 1995; Avery Sz Chiao 1996; 
Bachiller & Gutierrez 1997; Bachiller et al. 1998; Bachiller et al. 2001). In addition to the CH 3 OH 
lines, the C 34 S (J=2-l) line, the enhancement of which is also observed in some sources with 
outflow-core interaction (Sandell et al. 1994; Blake et al. 1995; Garay et al. 1998; Bachiller et al. 2001), 
was observed simultaneously. Belloche et al. (2002) also observed C 34 S ( J=2-l) as well as the main 
isotope in IRAM04191. 

2. Observations 

We carried out mapping observations of the CH3OH (Jk=2k-1r) and C S (J=2-l) lines 
with the 45 m telescope at Nobeyama Radio Observatory 3 from the 8th to the 11th of January, 
2000. The mapping observations have been performed at a grid spacing of 20" and a position 
angle of 30°, to align the mapping grid to the axis of the outflow. The observed molecular lines 
are summarized in Table 1. We used two SIS receivers with 500 MHz bandwidth by dividing the 
beam into two linear polarized waves through a wire grid, each of which covered the frequency 
range of the molecular lines listed in Table 1, and co- added the data to improve the signal to noise 
ratio. The beam squint between the two receivers was confirmed to be less than 1". Typical single 
sideband (SSB) system noise temperature was ~ 300 K for both receivers. The beam size and 
the main beam efficiency of the telescope was 17" and 0.5, respectively, at 96 GHz. The telescope 
pointing was checked every 60 - 90 minutes by observing SiO maser emission from NML-Tau at 43 
GHz, and was confirmed to be better than 5". The observational mode was position-switching. 

We used a bank of eight high-resolution acousto-optical spectrometers (AOS) with a bandwidth 
of 40 MHz and a frequency resolution of 37 kHz. The frequency resolution corresponds to ~ 0.11 
km s _1 for the observed lines. Typical rms noise per channel was ~ 0.06 K in units of Ta* (hereafter 
we present line intensities in units of Ta*). 

In order to calibrate the intensity variation caused by different receiver gains and/or sky 
conditions, we observed the central position of IRAM 04191 (02000 = 04 h 21 m 56 s . 9, (52ooo = 15° 
29' 46". 4) in every scan. We confirmed that the intensity variations during the whole observing 
period were less than 10 %. 



3 Nobeyama Radio Observatory (NRO) is a branch of the National Astronomical Observatory, an inter-university 
research institute operated by the Ministry of Education, Culture, Sports, Science and Technology of Japan. 
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3. Results 

3.1. CH 3 OH and C 34 S Spectra toward IRAM 04191 

Figure 1 shows profiles of the CH3OH and C 34 S lines toward the central position of IRAM 
04191. Since the integration time toward the central position was much longer (37 min) than that 
in other mapping positions, signal-to-noise ratios of these spectra (rms ~ 0.025 K) are better than 
those at other points. Three transitions of CH3OH, i.e., 2o-lo A + , 2_i-l_i E, and 2o~lo E were 
detected at the position of the protostar, while the 2!-li E transition was not detectable at any- 
observed positions, including the stellar position, at 1 a noise level of ~ 0.06 K. The C 34 S line was 
clearly detected toward the protostellar position, with intensity less than half of that of the CH3OH 
(2 -l A+) line. 

Both the CH3OH and C 34 S lines at the position of IRAM 04191 show single-peaked, Gaussian- 
like profiles with possible redshifted wings. Gaussian fittings to these profiles show that peak 
velocities and line widths of the spectra are ~ 6.6 - 6.7 km s _1 and ~ 0.4 - 0.5 km s _1 , respectively 
(see Table 2). The C 34 S profile is essentially same as in Belloche et al. (2002). Andre et al. (1999a) 
have shown that CS (J =2-1) and H 2 CO (2 12 - In) lines from IRAM 04191 exhibit the double- 
peaked, asymmetric spectral signature of protostellar collapse (Zhou 1992; Myers et al. 1995). The 
peak velocity of the CH3OH and C 34 S lines corresponds to the velocity of the absorption "dip" 
of the CS {J =2-1) and H 2 CO (2i 2 -ln) profiles, suggesting that the CH 3 OH and C 34 S lines are 
optically thin. Hereafter, we refer to 6.7 km s™ 1 as a systemic velocity of IRAM 04191, which is 
shown in a vertical dotted line in Figure 1. 

3.2. Spatial and Velocity Distribution of the CH3OH Emission 

In this subsection, spatial and velocity distribution of the CH3OH emission is discussed. We 
only adopt the CH3OH (2o-lo A + ) line here because the 2_i-l_i E transition shows essentially the 
same distribution as that the 2o-lo A + transition shows in our mapping observations, and because 
the 2o~lo E transition, which has a higher energy level than the other two detected lines (see 
Table 1), is detected only at several positions in the mapping observations. Figure 2 shows velocity 
channel maps of the CH3OH (2o~ lo A + ) line with a velocity interval of 0.2 km s -1 . The emission 
was significantly detected within the LSR velocity range from 5.0 to 7.4 km s" 1 . At blueshifted 
velocities ranging from 5.1 to 6.1 km s _1 , there are two CH3OH emission features detected: one 
appears ~ 60" south to the central star and the other at ~ 90" southwest to the star. The peak 
positions of these two emission features do not change much in this velocity range. At velocities 
from 6.3 to 7.3 km s _1 , around the systemic velocity, another CH3OH emission feature was detected 
in the vicinity of IRAM 04191. The peak position of this emission feature shifts from the northwest 
(at 6.3, 6.5 and 6.7 km s _1 ) to the southeast (at 6.9 and 7.1 km s" 1 ), and then to the northeast of 
the central star (at 7.3 km s^ 1 ). At redshifted velocities (6.9 - 7.3 km s" 1 ), there is an additional 
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CH3OH emission feature detected at ~ 60" west to the central star. This feature is likely to be 
independent of the 6.3 - 6.7 km s _1 component seen at the northwest, since the line profile toward 
this feature shows a double peak, one from the 6.3 - 6.7 km s _1 component and the other from this 
new feature (see Figure 3). At the northeastern edge of the mapping region, there appears strong 
CH3OH emission but the limited observing time prevented us from extending our mapping region 
to cover this component. Hence, this emission will not be discussed in this paper. We note that 
toward IRAS 04191+1523 (filled squares in Figure 2) the CH3OH emission distribution seems to 
be local minimum and there is no clear CH3OH condensation associated with IRAS 04191+1523. 

The above examination of the CH3OH channel maps shows that four CH3OH emission features 
at different velocities are present in IRAM 04191. In order to investigate the nature of these different 
emission features in more detail, the channel maps are integrated over three different velocity ranges, 
i.e., 5.0 - 6.5 km s _1 , 6.5 - 6.9 km s _1 , and 6.9 - 7.4 km s -1 , which roughly distinguish the different 
CH3OH emission components from each other. Figure 3 shows three maps integrated over the above 
velocity ranges, superposed on the CO (J=2-l) outflow map obtained by Andre et al. (1999a). 

The map of the systemic velocity (Figure 3; 6.5 - 6.9 km s" 1 ) shows a CH3OH condensation 
elongated in the east-west direction (P. A. ~ 110°), which is almost perpendicular to the outflow 
axis. The southern side of the condensation shows a concavity, through which the blueshifted 
outflow extends to south. This condensation is most likely a molecular envelope associated with 
IRAM 04191, although the position of IRAM 04191 is slightly (~ 20" east) offset from the cen- 
ter of the condensation. This kind of the elongated envelope is often associated with protostars 
(Mizuno et al. 1994; Zhou et al. 1994; Ohashi et al. 1997a; Ohashi et al. 1997b; Saito et al. 1999). 
Hence, hereafter we call this condensation "ENV" . The deconvolved size and the deconvolved line 
width of ENV were measured to be ~ 0.07 pc x 0.04 pc and ~ 0.5 km s _1 , respectively. These 
values are comparable to the typical size and line width of envelopes associated with low-mass 
protostellar sources (Benson Sz Myers 1989; Mizuno et al. 1994; Zhou et al. 1994). 

As has been suggested by the channel maps in Figure 2, there is a velocity gradient along the 
major axis of ENV. Figure 4, the position-velocity (P-V) diagram cutting along the major axis, 
demonstrates the velocity gradient clearly, as indicated by the thick dashed line. This velocity 
gradient is most probably attributed to rotation of ENV because it is obviously seen along the 
major axis of ENV. This velocity gradient was also seen in the N2H+ (J=l-0) line (Belloche et 
al. 2002). It should be noted that, the velocity gradient can be seen only in the vicinity of the 
central star (within ~ 0.01 pc from the central star), while the outer part of ENV, particularly 
the western side, does not show a clear velocity shift from the systemic velocity. This might be 
because the rotation of ENV is differential, so that the rotation velocity becomes larger in the inner 
part of ENV. The velocity gradient was measured to be ~ 40 km s _1 pc -1 , which is 10-100 times 
larger than that in other molecular cores measured by Goodman et al. (1993). This difference 
may result from the difference in the radius where velocity gradients are measured, i.e., in our case 
the gradient was measured at ~ 0.01 pc in radius, while Goodman et al. (1993) measured those 
gradients typically at ~ 0.1 pc in radius. 
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The map integrated over 5.0 - 6.5 km s^ 1 (the middle panel in Figure 3) shows two compact 
(~ 0.03 pc) condensations, connected by extended emission at lower levels. These condensations 
correspond to the two emission features seen at blueshifted velocities in the channel maps. Here- 
after, we will call these compact CH3OH condensations "BLUE1" (the western, stronger one) and 
"BLUE2", respectively. These two blueshifted gas components were also detected in the CS(2-1) 
line (Belloche et al. 2002; see their Fig. 6). As described below, physical properties of BLUE1 and 
BLUE2 are significantly different from those of ENV: 

(1) BLUE1 and BLUE2 are not located in the vicinity of the central star, but located at the 
southern edge of ENV. Comparison with the CO outflow map shows that BLUE1 and BLUE2 
locate on the the western and eastern sides of the blueshifted outflow lobe, respectively. 
Specifically, the peak of BLUE1 locates just beside the strongest peak in the blueshifted CO 
outflow. 

(2) Line profiles toward BLUE1 and BLUE2 show non-Gaussian shapes, with much wider line 
widths (up to 2.0 km s _1 ) than that of ENV and centroid velocities blueshifted by ~ 0.8 
km s _1 from the systemic velocity (Table 2, 3). The P-V diagram (Figure 5) of the CH3OH 
emission along the blueshifted outflow (solid line in Figure 3) shows abrupt increase of the 
line width from ~ 0.5 km s _1 to ~ 2.0 km s _1 at the position of BLUE 1. 

These results suggest that BLUE1 and BLUE2 are physically different from ENV. The origin of 
BLUE1 and BLUE2 will be discussed in §4 in detail. 

In the map integrated over 6.9 - 7.4 km s _1 (the lower panel in Figure 3), there is a CH3OH 
condensation detected to the west of the protostar, as well as a condensation appearing to the 
east of the protostar. As seen in Figure 2, the eastern condensation seems to be a part of ENV, 
while it delineates the southeast side of the redshifted outflow and could represent the interaction 
region between ENV and the redshifted outflow (Belloche et al. 2002). On the other hand, the 
western condensation is not a part of ENV, but the one corresponding to the emission feature 
detected at redshifted velocities in the channel maps. Hereafter, we will call the redshifted western 
condensation "RED". The line profile toward RED shows a double-peaked feature, i.e., one peaked 
at 6.7 km s _1 arising from ENV, and the other peaked at 7.1 km s -1 arising from RED. It is found 
that, while the size of RED is as compact as that of BLUE1 or BLUE2 (~ 0.03 pc), the line width 
is narrow (~ 0.3 km s^ 1 ). 

3.3. Spatial and Velocity Distribution of the C 34 S Emission 

Figure 6 shows velocity channel maps of the C S (J=2— 1) emission. The velocity range and 
interval of Figure 6 are same as those of the CH3OH (2o-lo A + ) channel maps shown in Figure 2. 
The typical intensity of the C 34 S emission is a factor of 2 lower than the intensity of the CH3OH 
emission in the mapped area. The C 34 S emission was not detected at every velocity channel where 
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the CH3OH emission was detected, but only in the velocity range from 6.3 to 7.1 km s where ENV 
and a part of RED were detected in the CH3OH map. In that velocity range, the spatial distribution 
of the C 34 S emission is similar to that of the CH3OH emission. Since the C 34 S emission is weaker 
than the CH3OH emission, we will concentrate on the CH3OH data in the following discussions. 



3.4. CH3OH Abundance 

The molecular abundance is one of the most important parameters to describe chemical prop- 
erties of the four condensations, that is, ENV, BLUE1, BLUE2, and RED. We have estimated the 
molecular abundance of CH3OH for each condensation by dividing the CH3OH column density 
(Nch 3 ok) by the molecular hydrogen column density (Nu 2 ), both of which are measured at the 
central position of each condensation. In order to estimate Ncn 3 ou, we assume that the emission 
lines are optically thin, which is valid for the emission lines we observed (see Table 3). In addi- 
tion we assume Local Thermal Equilibrium (LTE) condition of the CH3OH transitions and the 
equal partition function for both A-state and E-state CH 3 OH (Friberg et al. 1988). Under these 
conditions, the column density is given by the following formula; 

N = 8^^ Z(T ex ) T^AV 

CH3 ° H C 3 (2J, + 3)Aex P (-^)(l-exp(-^))J(T ex )-J(T bg )' U 

where 

J ( T ) = jt\ r ( 2 ) 

exp(^) - 1 

In the above expressions, h is the Planck constant, k is the Boltzmann constant, c is the speed of 
light, v is the line frequency, T ex is the excitation temperature of the CH3OH transitions, Tb g is the 
background radiation temperature, A is the Einstein A coefficient of the transition, Z is the partition 
function, E\ is the rotational energy level of the lower energy state, J\ is the rotational quantum 
number of the lower energy state, and Tg AV is the integrated line intensity. We adopted T cx of 10 K 
for all the condensations, according to our estimation of T ex from excitation analyses using the data 
of the (</k=5k~4k) band obtained with the JCMT and the data of the (Jk=2k-1k) band presented 
here on the assumption of the Large Velocity Gradient (LVG) condition (Goldreich & Kwan 1974; 
Scoville & Solomon 1974). Details of the LVG analyses with the multi-transitional CH3OH lines 
will be presented in a subsequent paper. 

The column densities of molecular hydrogen, on the other hand, were estimated using the 1.3 
mm dust continuum data of IRAM 04191 taken with the IRAM 30 m telescope (Andre et al. 1999b; 
Motte & Andre 2001). We read out intensity values of the dust continuum emission at each central 
position, and calculated iVH 2 by using the formula and the dust temperature used in Andre et al. 
(1999a) and Motte & Andre (2001). CH3OH abundances were estimated from these A^ch 3 oh and 
A^h 2 • Here, we did not take into account the difference of the beam size between our observations 
(17") and the 1.3 mm continuum observations with the IRAM 30 m telescope (13"). All the derived 
values are listed in Table 3. 
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The CH3OH abundance toward the protostar was estimated to be ~ 2.1 x 10 , which is 
indistinguishable from the CH3OH abundance toward the cyanopolyyne peak in TMC-1 (~ 2.0 x 
10~ 9 ; Pratap et al. 1997, Takakuwa et al. 2000). On the other hand, the CH3OH abundances 
toward BLUE1, BLUE2, and RED are 5-10 times larger than that toward the protostar. We will 
discuss this abundance enhancement later in more detail. 



3.5. Physical Properties of the CH3OH Condensations 

In addition to the chemical properties of the CH 3 OH condensations discussed above, we ex- 
amined the physical properties of the condensations. The LTE mass of each condensation, Mlte, 
was calculated using the column density and size listed in Table 3. The fractional abundance of 
CH3OH for each condensation we estimated in the previous section (also listed in Table 3) was 
used to estimate the total H2 gas mass. The LTE mass obtained in this way is 2.3 M for ENV, 
while less than 0.1 M for the other three condensations. In order to examine whether or not the 
CH3OH condensations are gravitationally bound, we compared Mlte with the virial mass (M vir ) 
calculated from the size (D) and the FWHM line width (AV) using the following formula, 

Mmr = (3) 



where 



AV 2 



C 2 ff ~ . (4) 

e// 81n2 K ' 

In the above expressions C e ff is the effective sound speed of gas and G is the gravitational constant. 

Note that the above expressions assume a condensation with a uniform density distribution. The 

masses of each CH3OH condensation derived from these two different ways are listed in Table 3. 

The virial mass of ENV is similar to its LTE mass, suggesting that ENV is probably virialized. 
In addition, the derived mass is comparable to the typical mass of the low-mass star-forming en- 
velopes (Mizuno et al. 1994; Zhou et al. 1994). These results suggest that ENV is an envelope being 
in the course of star- formation, and ENV is indeed associated with the protostar of IRAM 04191. 

In contrast to ENV, the virial masses of BLUE1 and BLUE2 are ~ 30-100 times larger than 
their LTE masses, suggesting that the blueshifted components are most likely gravitationally un- 
bound. Similarly to the blueshifted components, the virial mass of RED is ~10 times larger than 
its LTE mass, which suggests that RED is also likely gravitationally unbound. 
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4. Discussion 

4.1. The Origin of the Blueshifted and Redshifted CH3OH Condensations 

Our CH3OH observations of IRAM 04191 have revealed that, in addition to the envelope 
(ENV) associated with star formation activity, there are two blueshifted gas condensations (BLUE1, 
BLUE2) and one redshifted gas condensation (RED). The observed properties of the blueshifted 
components, i.e., located at the periphery of the blueshifted CO outflow, the wider line width, and 
the centroid velocity blueshifted from the systemic velocity, suggest that the blueshifted components 
are likely to be related to the blue lobe of the CO outflow. The total momentum of the blueshifted 
components is comparable to the CO outflow momentum (~ 0.12 km s _1 ; Andre et al. 1999a), 
suggesting that the blueshifted components can be pushed by the outflow. In addition, offset 
velocities (Voff = (Vch 3 oh ~~ Vsys) / cos i; where i is the inclination angle of the outflow from the 
line of sight) of the blueshifted CH3OH condensations are comparable to velocities necessary for 
the blueshifted condensations to escape from the gravitational potential of the central star-forming 
system, that is, V 9 = \J^§^-, where M is the mass of the envelope (i.e., = 2.3 M ) + star (0.05 
M Q ) and R is the projected distance from the protostar (Table 4). These results suggest that both 
B1UE1 and BLUE2 are most likely to be formed in consequence of interaction between the outflow 
and the ambient dense gas surrounding the protostar, and that the blueshifted components are 
pushed away from the ambient cloud core. The eastern condensation, which seems to be a part of 
ENV and delineates the southeast side of the redshifted outflow (Figure 3 lower panel), might also 
represent a core interacting with the redshifted outflow (Belloche et al. 2002). 

Nakano, Hasegawa, & Norman (1995) theoretically investigated mechanisms of setting the final 
stellar mass in the process of star formation, and proposed that the dispersal of the parent cloud 
core as a result of the interaction with the mass outflow (e.g. Momose et al. 1996; Velusamy k, 
Langer 1998) is crucial for fixing the final stellar mass. Hence, it is interesting to discuss whether the 
interaction with the IRAM 04191 molecular outflow plays an important role in the dissipation of the 
surrounding molecular cloud core and the process of the stellar mass determination. As discussed 
in §3.5, the blueshifted components are gravitationally unbound, suggesting that they will most 
likely dissipate eventually. The time scale of the dissipation, roughly estimated as D/C e ff where 
D is the size of the condensation and C e f f is the effective sound speed derived from the line width, 
is ~ 5 x 10 4 yr. This time scale is comparable to the typical life time of the accreting protostar, 
~ 10 5 yr (Adams et al. 1987; Adams 1990; Lada 1991). Hence, the dissipation of the surrounding 
cloud core may sufficiently proceed around IRAM 04191 during the main accretion phase, resulting 
in setting the stellar mass. 

We note, however, that the total mass of the blueshifted components (~ 0.15 M Q ), which will 
most likely dissipate, is at most 7% of the total mass of the envelope, suggesting that the total 
amount of the dissipating gas may not be high enough to have a significant impact onto the supply of 
matter for the mass accretion. Since IRAM 04191 is considered to be a very young class protostar 
(Andre et al. 1999a; Andre et al. 1999b), the outflow may just have started interacting with the 
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surrounding cloud core so that the core dissipation may not proceed enough yet. We speculate 
that the outflow will keep interacting with the surrounding cloud core, resulting in dissipation of 
significant fraction of the surrounding cloud core. 

The origin of the redshifted component, on the other hand, is less clear. Its narrow line width 
(~ 0.3 km s _1 ) suggests that the redshifted component is not related to the outflow. There is no 
clear star formation activity found in the redshifted condensation. The physical properties of the 
redshifted condensation, i.e., the narrow line width, the small mass (~ 0.07 M ) and size (~ 0.02 
pc), and non-virialization, are rather similar to those of small-scale clumps found in starless cores 
(Langer et al. 1995; Peng et al. 1998; Ohashi et al. 1999). 

4.2. CH3OH Abundance Enhancement 

Previous observations of molecular outflows associated with L1157 (Bachiller et al. 1995; Avery & Chiao 1996; 
Bachiller et al. 2001), NGC1333 IRAS 2 (Sandell et al. 1994; Bachiller et al. 1998), and BHR 71 
(Garay et al. 1998) have also showed CH3OH emission arising from regions where the outflows in- 
teract with the surrounding dense gas. An important, common result obtained from these previous 
observations is that the CH3OH emission arising from the interaction regions is much more intense 
and broader than that from the positions of the driving sources. The more intense emission has 
been attributed to bowshock heating, caused by the interaction between the outflows and the dense 
gas, leading to the enhancement of the CH3OH abundance (see Table 5). Is the blueshifted CH3OH 
emission in IRAM 04191 under influence of such bowshock heating? 

As discussed in §3.4., the CH3OH abundance in the blueshifted condensations is enhanced by 
a factor of ~ 10 with respect to that at the position of the protostar. This abundance enhancement 
might be due to shocks associated with the interaction between the blueshifted outflow and the 
surrounding dense gas. However, the enhancement seen in the blueshifted condensations is not as 
dramatic as that seen in the above-mentioned three outflows associated with shocks: the CH3OH 
abundances in these outflows with shocks were estimated to be enhanced by more than 2 orders 
of magnitude (Sandell et al. 1994; Bachiller et al. 1995; Bachiller et al. 1998; Bachiller et al. 2001; 
Avery Sz Chiao 1996; Garay et al. 1998). We also note that the C 34 S emission, which is also en- 
hanced in the shocked regions (Sandell et al. 1994; Blake et al. 1995; Garay et al. 1998; Bachiller et al. 2001) 
was barely detected at the positions of the blueshifted components in IRAM 04191. 

If CH3OH is enhanced due to shocks caused by the outflow, then the CH3OH enhancement 
should be localized on the shocked regions. In order to see if there is such localization of the CH3OH 
enhancement, the CH3OH abundance map in the IRAM 04191 region, shown in Figure 7, was made 
from the ratio of the total integrated intensity map of the CH3OH emission to the 1.3 mm dust 
continuum map (Andre et al. 1999a; Motte &; Andre 2001). We adopted the same assumptions as 
those in §3.4 to create this abundance map. We can clearly see that the CH3OH abundance is 
extensively enhanced in the periphery of the entire region, including the position of RED, which 
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shows no sign of interaction. There is no local enhancement of the CH3OH abundance toward 
BLUE2 and probably BLUE1, for which the fractional abundance was not determined because 
of less signal-to- noise ratio of the 1.3 mm dust continuum map (Motte & Andre 2001). We have 
also conducted sub-millimeter CH3OH (Jk=5k~4k, 7k~6k) observations with the JCMT toward 
BLUE1, and have found the temperature of the blueshifted component is cold (~ 10 K), suggesting 
that there is no local shock heating toward the blueshifted component (our forthcoming paper). 
These results suggest that the CH3OH abundance enhancement seen in the blueshifted components 
is not due to shocks associated with the interaction between the blueshifted outflow and the dense 
core. 

We consider that the observed abundance gradient is likely to reflect the chemistry in the 
protostellar core rather than the effect of shocks due to the outflow. Recently, Aikawa et al. (2001) 
have constructed chemical evolutionary models of a collapsing dense core. They found that chemical 
abundances of molecules are functions of both the evolutionary time-scale and the distance from 
the core center. They showed that, the abundances of "early-type" molecules, which are considered 
to be more abundant in the early stage of chemical evolution of dense cores, significantly decrease 
in the central regions of the dense core as the evolutionary time goes on, resulting in abundances 
of early-type molecules appearing "enhanced" in the periphery of the dense core. Since CH3OH 
is considered to be one of the early-type molecules (Pratap et al. 1997; Takakuwa et al. 1998; 
Takakuwa et al. 2000), the observed CH3OH abundance variation in IRAM 04191 is consistent 
with their model. Such abundance gradients of early-type molecules such as CCS were also shown 
in other observations (Kuiper et al. 1996; Ohashi et al. 1999; Ohashi 2000). 

It is interesting to discuss why the blueshifted CH3OH emission in IRAM 04191 is not due to 
shock heating by the blueshifted outflow, which most probably interacts with the ambient cloud. 
In Table 5, dynamical properties of the CO outflow from IRAM 04191 are compared with those of 
the three outflows with shock-enhanced CH3OH emission. The table shows that the flow velocities 
(Vfl ow ) of the CO outflows are very different between the three CO outflow and IRAM 04191 outflow. 
The flow velocities of the three CO outflows in the shocked regions are all above 20 km s _1 , whereas 
the flow velocity of IRAM 04191 outflow is below 10 km s -1 . We also note that the line widths and 
the offset velocity (V ff) of the CH3OH components in IRAM 04191 are also small as compared 
to those in the other shocked outflows (Table 5). Shocks associated with molecular outflows from 
protostars are generally considered to be C-type shocks (Hollenbach 1997). According to detailed 
models of C-type shocks including the calculation of the chemical composition (i.e., coolants) and 
the grain size distribution, heating caused by C-type shocks is a strong function of the flow velocity, 
or the shock velocity (Kaufman & Neufeld 1996a; Kaufman & Neufeld 1996b; Bergin et al. 1998). 
These models show that shock velocity higher than ~ 10 km s _1 is necessary to yield temperature 
enhancement higher than 120 K, which is considered to be the minimum required temperature to 
release CH3OH on dust grains into gas phase (Sandford & Allamandola 1993). Thus, the difference 
in Vflow may explain why IRAM 04191 is not associated with CH3OH emission enhanced due to 
shocks. 
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5. Conclusions 

We have carried out mapping observations around the young Class protostar, IRAM 04191+1522, 
in the CH 3 OH ( J k =2 k -1k) and C 34 S ( J =2-1) lines with the NRO 45 m telescope. The main results 
are as follows. 

1. A dense gas envelope surrounding the protostar was traced both in the CH 3 OH and C 34 S 
lines. The envelope has a size of ~ 0.07 pc x 0.04 pc, line width of ~ 0.5 km s^ 1 , and LTE mass 
of ~ 2.3 M Q . There is a steep velocity gradient of ~ 40 km s _1 pc -1 along the NW-SE direction, 
i.e., perpendicular to the axis of the outflow, in the envelope. This velocity gradient is interpreted 
as the rotation of the envelope. 

2. The CH3OH map shows two compact (~ 0.03 pc in size), blueshifted (~ 0.8 km s _1 from 
the systemic velocity) condensations with wide line widths (up to 2.0 km s _1 ) at the periphery of 
the blueshifted lobe of the CO outflow. The compact condensations have momenta (~ 0.06 M 
km s _1 ) comparable to that of the blueshifted molecular outflow. These results suggest that the 
two blueshifted components are likely to be formed by means of interaction between the outflow 
and its surrounding dense gas. The LTE masses of these blueshifted components (~ 0.07 M ) are 
two orders of magnitude lower than their virial masses, suggesting that these condensations are not 
bound by their self-gravity, and most likely will dissipate eventually. 

3. The CH3OH abundance in the two blueshifted condensations is enhanced by a factor of ~ 10 
with respect to that toward the position of the protostar (~ 2.1 x 10~ 9 ). However, this abundance 
enhancement is less significant as compared to that seen in the outflows with shock heating. In 
addition, the spatial variation of the CH3OH abundance does not show the local enhancement at the 
positions of the blueshifted condensations, but shows the trend of decrease toward the center and 
increase toward the periphery. Such an abundance variation is explained by the effect of chemical 
evolution of the collapsing dense core rather than by the shock chemistry. 

4. A comparison of the properties of the CO outflow in IRAM 04191 with those with shock 
heating suggests that the velocity of the outflow in IRAM 04191 (less than 10 km s _1 ) is not high 
enough to produce shocks that can release the CH3OH molecules from dust mantles. 

We are grateful to M. Choi for his fruitful discussions. We would like to thank all the NRO staff 
supporting this work. S.T. was supported by a postdoctoral fellowship of the Institute of Astronomy 
and Astrophysics, Academia Sinica. N.O. is supported in part by NSC grant 91-2112-M-001-029. 
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Fig. 1.— CH3OH (J K =2 -1 A+, 2_i-l_i E, 2o-l E, and 2i-li E) and C 34 S (J=2-l) line 
profiles toward the protostellar position of IRAM 04191+1522. The vertical dotted line indicates 
the systemic velocity of 6.7 km s -1 . 

Fig. 2.— Velocity channel maps of the CH 3 OH (J K =2 -1 A + ) emission in IRAM 04191+1522. 
The channel widths are 0.2 km s^ 1 and the center velocities are shown in the Figure. Contour 
levels are from 0.12 K by steps of 0.08 K. Cross marks indicate the position of the protostar, IRAM 
04191+1522. Filled squares shows the position of IRAS 04191+1523. Dot marks are the observed 
points. 

Fig. 3. — Velocity channel maps of the CH3OH ( Jk=2o-1o A + ) emission in three velocity regimes 
(Systemic : 6.5 - 6.9 km s -1 , BLUE : 5.0 - 6.5 km s _1 , and RED : 6.9 - 7.4 km s _1 ) plus line profiles 
toward the representative positions (ENV, BLUE1, BLUE2, and RED) in IRAM 04191+1522. 
Contour levels in the channel maps are from 0.10 K km s _1 by steps of 0.04 K km s -1 . In the 
channel maps we overlay the map of the CO outflow by Andre et al. (1999a) (Red contours indicate 
red lobe and blue contours blue lobe of the CO outflow). Marks in the channel maps have the same 
meanings as in Figure 2. Straight lines in the channel maps indicate the cut of the position-velocity 
diagram shown in Figure 5. In the line profiles we show the systemic velocity of 6.7 km s _1 by a 
dashed vertical line. 

Fig. 4. — A position-velocity (P-V) diagram of the CH3OH (Jk=2o-1o A + ) emission along the 
major axis of the protostellar envelope (P.A.=120°) in IRAM 04191+1522. Contour levels are from 
0.12 K by steps of 0.08 K. A vertical dotted line in the figure indicates the systemic velocity (= 6.7 
km s _1 ), and a horizontal line in the figure represents the position of the protostar. A bold dotted 
line indicates the velocity gradient of the protostellar envelope around the protostar. 

Fig. 5. — A P-V diagram of the CH3OH (Jk=2o-1o A + ) emission along the cut which crosses 
BLUE1 with P.A.=30° (straight lines in Figure 3). Contour levels are from 0.15 K by steps of 0.10 
K. We can see sudden increase of the line width at the position of BLUE 1. 

Fig. 6. — Same as Figure 2 but for the C 34 S ( J =2-1) emission. 

Fig. 7. — A map of the CH3OH abundance derived from the ratio of the total integrated intensity 
map of the CH3OH (Jk=2o~1o A + ) line to the 1.3 mm dust continuum map by Motte & Andre 
(2001) in IRAM 04191+1522. The innermost contour corresponds to the CH3OH abundance of 
1.6 x 10 -9 , and the abundance increases by a factor of 2, 3, 4, 10 toward the outer contours as 
indicated in the figure. In the figure we plot the position of IRAS 04191+1523 (north-eastern), 
IRAM 04191+1522 (central), BLUE1 (south-western ; outside the border of the map), BLUE2 
(south of IRAM 04191+1522), and RED (west of IRAM 04191+1522) by cross marks. The signal- 
to-noise ratio of the 1.3 mm dust continuum map toward BLUE1 is worse than the other positions 
and therefore BLUE1 is outside the map boundary. We can see that the abundance systematically 
increases in the outer part of the protostellar envelope. 
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Table 1. Observed Molecular Lines 



Molecule 


Transition 


Frequency 
(GHz) 


Einstein A Coefficient 
(s- 1 ) 


E« 
(K) 


CH 3 OH a 


J K = 2 -l A+ 


96.74142 


3.41 x 10~ 6 


6.96 c 




Jk = 2_i-l_i E 


96.73939 


2.56 x 10 -6 


4.64 c 




Jk = 2o-lo E 


96.74458 


3.41 x 1(T 6 


12.2 C 




J K = 2i-li E 


96.75551 


2.62 x 10" 6 


20.1 c 


C 34 S b 


J = 2-1 


96.412982 


1.60 x 10" 5 


6.93 



a Line parameters of CH3OH are from Xu & Lovas (1997). 

b Line parameters of C 34 S are from Swade (1989). 

c Energy level from the ground state of the relevant symmetric state. 



Table 2. Observed Line Parameters of CH 3 OH and C 34 S 



CH 3 OH a C 34 S b 







Jk = 2o 


-lo A+ 




J K = 2-1-1. 


-1 E 




Jk = 2o- 


lo E 




J = 2-1 




ition 


T A * 


AV 


"V peak 


T A * 


AV 


'V peak 


T A * 


AV 


'V peak 


T A * 


AV 


Vpeafc 




(K) 


(km s _1 


) (km s- 1 ) 


(K) 


(km s _1 ) 


(km s _1 ) 


(K) 


(km s _1 ) 


(km s _1 ) 


(K) 


(km s _1 ) 


(km s _1 ) 


NV C 


1.207 


0.54 


6.70 


0.909 


0.53 


6.63 


0.225 


0.43 


6.65 


0.446 


0.50 


6.66 


UEl d 


0.599 


0.85 


6.12 


0.519 


1.02 


5.98 


0.272 


0.21 


5.87 


0.185 


0.36 


5.26 


UE2 d 


0.283 


1.61 


6.50 


0.234 


1.96 


5.48 














ED C 


0.740 


0.38 


6.72 


0.544 


0.48 


6.69 








0.208 


0.24 


6.55 


ED C 


0.908 


0.37 


7.14 


0.639 


0.32 


7.07 








0.245 


0.45 


6.95 



ms noise levels of the CH 3 OH lines are 0.024 K, 0.059 K, 0.041 K, and 0.052 K in ENV, BLUE1, BLUE2, and RED, respectively. 
,ms noise levels of the C 34 S line are 0.026 K, 0.058 K, 0.055 K, and 0.062 K in ENV, BLUE1, BLUE2, and RED, respectively, 
ine parameters are determined by a single Gaussian fit to the spectra. 

ine parameters are determined by measuring the peak intensity, peak velocity, and the FWHM of the spectra, 
ine parameters are determined by two-components Gaussian fit to the spectra. 



Table 3. Physical and Chemical Properties of the CH3OH Components 



nponent 




T b,c 


^CH 3 OH C ' d 


V c,d 
-*-CH 3 OH 


VcH 3 OH e 


#FWHM f 


D { 


AV C 




M LTE g 


Mom.S' h 




(x 10 22 cm" 2 ) 




(x 10 14 cm' 2 ) 


(x 10~ 8 ) 


(km s _1 ) 


(" X ") 


(PC) 


(km s _1 ) 


(Me) 


(M ) 


(M km s" 1 ) 


ENV 


4.26 


0.43 


0.89 


0.2 


6.7 


106 x 51 


0.050 


0.53 


2.4 


2.300 




iLUEl 


0.54 


0.19 


1.16 


2.1 


5.9 


27 x 21 


0.016 


0.84 


1.5 


0.064 


0.058 


SLUE2 


0.57 


0.09 


0.62 


1.1 


6.0 


74 x 26 


0.030 


1.61 


8.6 


0.085 


0.073 


RED 


0.73 


0.31 


0.83 


1.1 


7.1 


56 x 19 


0.022 


0.35 


0.7 


0.069 


0.040 



lolecular hydrogen column densities are estimated from the 1.3 mm dust continuum map by Motte & Andre (2001). 

cak optical depth of the 2o-lo A + line. 

)n the assumption of the excitation temperature of 10 K. 

'n the assumption of equal abundance of A- and E-symmetric species of CH3OH. 

picked up the CH3OH spectra toward the central positions of each identified component, and made Gaussian fittings to the spectra to measure 
'elocity and the line width (AV) of each gas component. If the line profiles are not well-fitted by a Gaussian, we simply measure the FWHM 
width. In estimating the line widths, we deconvolved the measured line width with the velocity resolution. 

T e measured the size of each gas components by 2-dimensional Gaussian fittings to the integrated intensity maps in Figure 3. The deconvolution 
the beam size were performed to obtain the size. Here, the size (D) of the condensations is the geometrical mean of the major and minor axis 
hm) from the 2-dimensional Gaussian fit. 

3e text. 

orrections for the inclination of the outflow to the line of sight (~ 40° : Andre et al. 1999a) are applied. 
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Table 4. Offset Velocities of the CH 3 OH Condensations 



Condensations R V // V g 

(pc) (km s _1 ) (km s _1 ) 



BLUE1 


0.055 


1.03 


0.60 


BLUE2 


0.027 


0.87 


0.86 


RED 


0.041 


0.45 a 


0.71 



a No correction for the inclination angle is ap- 
plied. 
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Table 5. Comparison of Outflow Sources with CH3OH Observations 





Protostar 




CO Outflow 






CH3OH 






Source 


L&oi 




V flow a '° 


Momentum 13 


T dyn 


X c 


2 -lo A+ AV 


Vo// D 


Rcf. d 




(L©) 


(K) 


(km s- 1 ) 


(M km s- 1 ) 


(y) 




(km s- 1 ) 


(km s- 1 ) 




IRAM 04191 


0.15 


18 


5.5 


0.12 


8000 


<10 


0.85 


1.0 


1,2,3 


L1157 


5.8 


42 


49 


30 


15000 


400 


5.9 


36 


4,5,6 


NGC1333 IRAS 2 


7.2 


26 


22 


0.64 


12000 


300 


7.2 


18 


7,8,9 


BHR 71 


9 


35 


28 


11 


10000 


40 


4.1 


31 


10,11 



a The flow velocity (Vfl ow ) is determined by dividing the momentum by the outflow mass, that is, the flow velocity weighted by the 
mass. 

b Corrections for the effect of the inclination of the outflow are applied. 

C CH30H abundance enhancement toward the outflow as compared to that in the protostcllar position. 

d References : (1) This work, (2) Andre ct al. 1999a, (3) Ohashi 1991, (4) Bachiller ct al. 2001, (5) Shirley et al. 2000, (6) Bachiller 
ct al. 1995, (7) Knee & Sandell 2000, (8) Bachiller et al. 1998, (9) Sandcll ct al. 1994, (10) Garay ct al. 1998, (11) Bourke et al. 1997 
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